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^ ■ ABSTRACT 

(/"■) . The conversion of gas into stars is a fundamental process in astrophysics and cosmol- 

£N| ■ ogy. Stars are known to form from the gravitational collapse of dense clumps in interstellar 

molecular clouds, and it has been proposed that the resulting star formation rate is propor- 
qq . tional to either the amount of mass above a threshold gas surface density, or the gas volume 

1 density. These star-formation prescriptions appear to hold in nearby molecular clouds in our 

fvq \ Milky Way Galaxy's disk as well as in distant galaxies where the star formation rates are of- 

ten much larger. The inner 500 pc of our Galaxy, the Central Molecular Zone (CMZ), contains 
the largest concentration of dense, high-surface density molecular gas in the Milky Way, pro- 
viding an environment where the validity of star-formation prescriptions can be tested. Here 
■ we show that by several measures, the current star formation rate in the CMZ is an order-of- 

$_i \ magnitude lower than the rates predicted by the currently accepted prescriptions. In particular, 

the region 1° < I < 3.5°, \b\ < 0.5° contains ~ 10 7 M Q of dense molecular gas — enough to 
form 1000 Orion-like clusters — but the present-day star formation rate within this gas is only 
equivalent to that in Orion. In addition to density, another property of molecular clouds, such 
as the amplitude of turbulent motions, must be included in the star-formation prescription to 
predict the star formation rate in a given mass of molecular gas. 

Key words: stars formation, ISM:evolution, radio lines:ISM, line:profiles, masers, stars:early 
type 



1 INTRODUCTION 

Stars play a pivotal role in shaping the cosmos. High-mass stars 
contributed to the ionisation of the Universe after the Cosmic Dark 
Ages. They drive energy cycles and chemical enrichment on galac- 
tic scales, hence sculpting galactic structure. Thus, the conversion 
of gas into stars is fundamental to astrophysics and cosmology. 
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The rate at which gas is converted into stars has been mea- 
sured in the disks of nearby galaxies. When averaged over hun- 
dreds of parsecs, the star-formation rate (SFR) was found to have 
a power-law dependence on the gas sur face-density as described 
by the Schmidt-Kennicutt (SK) relations JSchmidtl 1 959l : iKennicuttl 
ll998l : lKennicutt & Evansll2012h . A linear relationship is found be- 
tween SFR and gas surface-density above a local extinction thresh- 
old, Ak ~0.8 magnitudes at a near-infrared wavelength of 2.2 
/xm, corresponding to a gas column-density of ~ 7.4 x 10 21 
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hydrogen molecu les per cm 2 or a gas surface-density E gas , of 
~ 116M©pc~ 2 dLada etalj|2012h . For a typical 0.1 pc radius 
cloud core, this corresponds to a volume density 3 x 10 4 hydrogen 
molecules cm -3 . 

Alternatively, it has been proposed that the SFR density is lin- 
early proportional to the mean gas density divided by the free-fall 
time multiplied by an efficiency factor estimated to be about 1% 
jKrumholzetai]|2012h . These surface- and volume-density rela- 
tions potentially unify our understanding of SFRs from the nearest 
star forming regions to ultra-luminous infrared galaxies (ULIRGs) 
and even star forming regions at intermed iate to high redshifts 
dSwinbank et alj20ld : lDanielson et alj|201 ll) . 

With the aim of testing these star formation (SF) scaling re- 
lations across large volumes of the Milky Way (MW), in § [2] we 
start by looking at the choice of observational diagnostics of dense 
gas and star formation activity in the Galaxy. In §[3] we present the 
survey data used in this work. In § [4] &|5]we discuss the striking 
differences between the dense gas and SF activity tracers across 
the Galaxy and investigate if any observational or systematic bi- 
ases may be causing such differences. We find no evidence that the 
results are strongly affected by observational or systematic biases. 
In §[6] we then perform a quantitative analysis of the SFR and gas 
mass in the CMZ and directly compare this to predictions of pro- 
posed column density threshold and volumetric SF relations. In §|7] 
we summarise the results and discuss the implications for the uni- 
versality of star formation relations. 



2 OBSERVATIONAL DIAGNOSTICS OF DENSE GAS 
AND STAR FORMATION ACTIVITY IN THE MILKY 

WAY 

We aim to test star formation scaling relations for representative 
samples of SF regions across large volumes of the MW. Such sys- 
tematic studies are now possible thanks to the (near) completion of 
many blind, multi-wavelength, continuum and spectral line Galac- 
tic plane surveys. However, care must be taken in the choice of 
SF activity diagnostics. All surveys aiming to observe a large frac- 
tion of the dense gas in the Galaxy must look directly through the 
Galactic disk plane. Extinction therefore makes traditional UV, op- 
tical and even near-IR diagnostics of SF activity unusable for all 
but the closest regions. 24/^m emission, commonly used to trace 
SF in extragalactic studies, suffers from extinction only in extreme 
environments (like the GC) but contamination from evolved stars 
can be significant, requiring other wavelength data to determine the 
source nature. Sensitivity and completeness limits mean counting 
young stellar objects still heavily embedded in their natal molecu- 
lar clouds can only be done reliably for the closest SF regions. 

Diagnostics are required at wavelengths longer than a few 
hundred microns, for which the Galaxy is effectively optically-thin. 
Sub-mm dust continuum emission is a robust gas mass tracer, but 
from this alone it is difficult to impose a volume density cutoff on 
the sample selection. Molecular line emission has the advantage 
that it can be used to probe gas within a certain density range. The 
actual density of gas probed by observations of a given transition is 
related to the critical density of that transition (the density at which 
the collisional de-excitation rate equals the Einstein A coefficient), 
but also depends non trivially on the sensitivity of the observations, 
the line opacity and gas kinetic temperature. A useful conce pt, tak- 
ing in to account these affects, is the effective critical density dEvans! 
1 19991) . For a given transition, this is the density of gas probed by ob- 
servations of a particular brightness temperature sensitivity for gas 



at a given kinetic temperature. Assuming detected emission comes 
from gas at a similar kinetic temperature (a reasonable assumption 
for dense gas clumps in molecular clouds) and that beam dilution is 
not an issue, a uniform sensitivity molecular line transition survey 
offers a natural way to probe the amount of gas at, and above, the 
effective critical density. 

Masers are well known indicators of SF activity and with 
strong transitions in the cm/mm wavelength regime they can be 
detected at large distances. However, it is not currently possible to 
derive an absolute star formation rate from first principles using 
masers. Free-free emission from ionised gas in Hn regions created 
by young high mass stars is readily traced by cm-continuum and 
radio recombination line (RRL) emission, both of whi ch are read- 
ily ob servable over large wavelength/frequency ranges dPeters et al.l 
l2012h . The number count of Hn regions can give a handle on the 
relative numbers of SF sites across the Galaxy. The absolute star 
formation rate can be estimated from the luminosity of the cm- 
continuum emission by calculating the rate o f ionising photons 
from high-mass stars dMurrav &~R ahman 2010), although contam- 
ination from non-thermal emission (which can lead to an overesti- 
mate of the SF activity) must be taken into account, especially to- 
wards the GC. It should be noted that each of these SF diagnostics 
traces the SF activity over different timescales. Masers are only ob- 
served towards regions that are still actively forming stars, so sign- 
post SF activity on timescales of a few 10 5 yr. Free-free emission 
from Hn regions is observed over the lifetimes of the high mass 
stars providing the ionising photons ie. timescales of a few Myr 
dMurrav & Rahmaij|201ol) . When we refer to the SFR throughout 
this work, we are implicitly talking about star formation activity 
over these timescales, which is significantly shorter than other di- 
agnostics (e.g. UV emission). 



Of the rec ent and planned survey; 
Survey (HOPS: IWalshe t al. 200& 1201 li - 



the H2O Galactic Plane 



Purcell et alj|2012h is the 



only blind molecular line survey covering a large fraction (100 
deg 2 ) of the Galactic plane simultaneously in multiple thermal 
molecular lines (including the important NH3 molecule), masers 
and RRLs . Below we describe HOPS and several newly available 
Galactic plane surveys (see Table [T} which we use to investigate 
how the SFR relations and proposed density thresholds hold for 
more representative SF regions across the MW. 



3 OBSERVATIONAL DATA 

Determining the rate of SF per unit mass of gas requires measuring, 
i) the total gas mass and, ii) how much SF is taking place within 
this gas. We make use of recently available NH3(1,1) data from 
HOPS and 70 to 500 /i m data from the Hers chel Infrared Galactic 
Plane Survey (HIGAL: iMolinari et ai1l20ld) to trace the mass of 
dense gas across the Galaxy. Ammonia (NH3) is well known as an 
excellent tracer of de nse molecular gas a t all evolutionary phases 
of the SF process dHo & Townesl 1 19831: ILongmore et al.l [20071 : 
iBeufher et al.l2009l : lHill et al.l20ld : lLongmore et al.l201ll) . includ- 
ing the earliest phases where low temperatures mean molecules like 
CO freeze out on to dust grains and no longer reliably trace the gas. 
Given the HOPS sensitivity of 0.2 K we expect to detect gas with 
an effective critical density of a few 10 3 , probin g gas close t o, and 
above, the proposed volume density threshold dEvans|[T999l) . The 



1 Although other large Galactic, molecular-line mapping proj e cts ar e 
also being conduc t ed on the Mopra telescope e.g. iBarnes et alj J201ll) ; 



iFoster et al] fepl ll) ; |jones et alj 420121) 
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500 /im HiGAL data traces optically-thin dust emission, provid- 
ing a measurement of the gas column density along a given line 
of sight. To directly compare the 500/^m and NH3(1,1) emission, 
we first re-binned the 35" 500/im map to the HOPS pixel scale 
using Montage (http://montage.ipac.caltech.edu/), then applied the 
HOPS NHs(l,l) emission mask to the re-binned HiGAL map. This 
provided the 500/^m emission from all HOPS NH3 emission re- 
gions. 

We use data from HO P S and the methano l multib eam (MMB 
ICaswell et all l20ld, 1201 ll : iGreen et alj I201GL 120121) surveys to 
detect the 22 GHz water maser emission and 6.7 GHz Class II 
methanol maser emission, respectively, across the Galactic plane. 
Both water and methanol masers are well known indicators of SF 
activity within molecular clouds. In particular, the 6.7 GHz Class 
II methanol maser trans ition is known to e xclusively trace regions 
of high mass (>8 M Q ) dMinier et al.ll200#l . Water masers are de- 
tected towards evolved stars, but these are found outside the disk- 
plane and the emissio n is typically two orders of magnitude weaker 
dCaswell et ail fcoiO). So sensitivity-limited, water maser Galactic 
plane surveys should predominantly trace SF regions. Removing 
evolved star contaminants is typically s traightforward using other 
readily-available multi-wavelength data dLumsden et alj|2002fl . 

We use the radio recombination line (RRL) data from HOPS 
and the Gre en Bank Telescope (G BT) H11 Region Discovery Sur- 
vey (HRDS: Anderson et al. 201 1) to compare the Galactic Hn re- 
gion number distribution to the dense gas and maser distributions. 
The HRDS completeness limit is sufficient to detect all Hn re- 
gions ionised by single O-stars to a distance of 12 kpc. Both HOPS 
and HRDS use RRLs to detect emission from Hn regions and the 
kinematic dimension helps to mitigate line-of-sight confusion that 
can affect continuum surveys. HRDS is much more sensitive than 
HOPS but is a targeted survey based on previous large-area cm- 
continuum and IR surveys. The HOPS RRL emission is less sen- 
sitive but allows for direct comparison with the maser, NH3 and 
500 /im data over the full 100 deg 2 of the HOPS survey region. 

The Hn region number counts provide a measure of the rela- 
tive distribution of SF sites acros s the Galaxy but we m ake use of 
recent analysis of WMAP data dMurrav & Rahmanll2010h to derive 
the absolute SFR in the GC. 



4 COMPARING DENSE GAS EMISSION AND SF 
ACTIVITY TRACERS ACROSS THE GALAXY 

As the Galaxy is optically-thin to the NH3(1,1), 500 /im, water 
maser, methanol maser and RRL emission, the emission surface 
density ratios between these tracers in the plane of the sky is equiv- 
alent to that in a face-on view of the Galaxy. Comparing the emis- 
sion ratios as a function of longitude and making the reasonable as- 
sumptions that i) the integrated intensity of the far-IR and NH3(1,1) 
emission is proportional to the mass of dense gas, and, ii) the num- 
ber of water masers, methanol masers, Hn regions and RRL in- 
tegrated intensity is proportional to the amount of SF, we now in- 
vestigate the relationship between the dense gas and SFR in the 
Galaxy. 

Figure Q] shows the dense gas emission and current star for- 
mation activity tracers as a function of Galactic longitude (I) and 
latitude (b). The dense gas distribution is dominated by the very 

2 Although the origin of the masers in disks or outflows is deba ted (e.g. 
lNorrisetal.h9 98; Minier et al. 2000; lDe Buizer et al.l2009ll2012h . 



bright and spatially-extende d emission within a few degrees of the 
Galac tic centre - the CMZ dMorris & Serabvn|[l996l ; lFerriere et al.1 
2007). This is easily distinguished from the rest of the Galactic 
plane by the very high intensity of dense gas emission and large gas 
velocity dispersion as indicated by wide spectral lines. While the 
dense gas emission is highly concentrated in the CMZ, the distri- 
bution of star formation activity tracers is relatively uniform across 
the Galaxy. Quantitatively, the CMZ accounts for ~80% of the inte- 
grated NH3(1,1) intensity but only contains 4% of the survey area. 
Yet, the CMZ does not stand out in either water or methanol maser 
emission, or in the number of Hn regions, which all trace recent 
hi gh-mass star formatio n. A qualitatively similar trend is reported 
bv lBeuther etalH2012h who compare the sub-mm dust continuum 
emission and GLIMPSE point sources as a function of Galactic 
longitude. 

Figure [2] shows Galactic longitude distributions of dense gas 
and star-formation activity indicators summed over the observed 
latitude range (—0.5° ^ b ^ 0.5°) for each longitude pixel. To 
make a direct comparison of dense gas and SF tracers as a function 
of longitude, we first sampled in 2-degree longitude bins to en- 
sure the volume of the Galaxy covered in each bin contains a large 
number of SF region s and so is appropriate for testing SF relations 
( IOnoderaet"aLll2 010). The emission or number count in each 2 de- 
gree longitude bin was then normalised by the total emission or 
number in the full survey area. 

The top panel of Figure [3] shows that there is a strong cor- 
relation between independent dense gas tracers. To separate the 
CMZ emission from that in the rest of the Milky Way, from here 
on we refer to the region \l\ < 5° with NH3(1,1) line widths 
AV > 15kms _1 as "GC-only", to distinguish it from the rest 
of the Galaxy, which we refer to as the "non-GC" region. No off- 
set is seen in the correlation between the GC-only and non-GC re- 
gions. The middle panel of Figure [3] shows a correlation between 
the dense gas and SF tracers for the non-GC regions. However, the 
GC-only regions are clearly distinct, with at least an order of mag- 
nitude brighter dense gas emission for the number of SF activity 
tracers. 

The bottom panel of Figure [3] shows the resulting dense gas 
vs SF tracer surface density ratio in Galacto-centric radius, Rgc, 
annuli of 0. 5 kpc. Rgc was calc ulated using the Galactic rota- 
tion curve of iBrand & Blitj d 19931) and assuming a distance to the 
Galactic centre of 8.5 kpc and a solar velocity of 220 kms~ . The 
gas between the CMZ and Rgc ~3 kpc shows emission at anoma- 
lous velocities so the rotation curve does not place reliable con- 
straints on Rgc over this range. The surface density ratios over this 
region, highlighted by the hatched rectangle, should be ignored. 
The ratio is approximately constant at Rgc >3 kpc. This suggests 
the linear relationship observed between the quan tity of gas above 
the proposed extinction threshold and the SFR dGao & Solomon! 
2004; Wu et al.ll2005ULada et al .11201 2) extends to a larger number 
of more representative SF regions across the Galaxy. By compari- 
son, the dense gas surface density towards the GC (RGC<0.5kpc) 
is orders of magnitude too large compared to the SF activity surface 
density. 



5 VARIATION IN SFR PER UNIT MASS OF DENSE GAS 
BETWEEN THE DISK AND GALACTIC CENTRE? 

We now investigate the striking difference between the longitude 
and surface density distribution of dense gas and SF tracers be- 
tween the GC-only and non-GC regions (Figure |3j- If the number 
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of masers and Hn regions were directly proportional to the inte- 
grated intensity of the dense gas emission, as would be expected 
if assumptions i) and ii) outlined in §[4]hold, and the SFR was set 
by the amount of dense gas, how many would be expected towards 
the GC-only region? With ~ 1/25 of the HOPS survey area and 4 
times the total NH3(1,1) integrated intensity outside the GC, the 
GC-only region has ~100 times the dense gas integrated intensity 
per unit survey area compared to the non-GC region. Therefore, the 
expected number of masers and Hn regions towards the GC should 
also be 100 times larger per unit area. The SF activity per unit area 
as traced by methanol masers is a factor 3 to 6 high er in the inner 
250 p c of the Galaxy compared to the spiral arms dCaswell et al.l 
|2010|) . However, our results show that this increase in maser surface 
density is dwarfed by the increase in the dense gas surface density. 
Given the amount of dense gas, the GC-only region appears defi- 
cient in SF tracers by two orders of magnitude. Taken at face value 
this suggests large differences in SFR per unit mass of dense gas 
between the disk and Galactic centre. Below we investigate if such 
a large deficit could be caused by observational or systematic bi- 
ases affecting the emission from one or more of the tracers used in 
this study. 

5.1 Potential observational biases 

These data are from blind large-area surveys with approximately 
uniform sensitivity. Therefore, the results should not be affected by 
target selection criteria or variable sensitivity. But what about the 
sensitivity limits themselves? For example, could HOPS be miss- 
ing a large fraction of the dense molecular gas outside the GC- 
only regio n as this NH3(1,1) e mission falls beneath the sensitivity 
threshold? Furcell et al.l Uovh show HOPS should detect 400 M 
clumps to 5.1 kpc and 5000 Mq clumps to the Galactic centre dis- 
tance. These limits suggest completeness is not a problem. 

Alternatively, as the GC-only clouds are more distant on av- 
erage than the non-GC regions, could there be a large number of 
weak water masers towards the GC-only region which fall below 
the HOPS detection limit? Ott et al (in prep.) have mapped the en- 
tire CMZ from 20 to 28 GHz with Mopra in a similar setup to the 
HOPS observations but with ~ 8 x better sensitivity and using the 
correlator in broadband rather than zoom mode (resulting in veloc- 
ity resolution of 3.5 kms -1 rather than 0.4 kms -1 as for HOPS). 
They find approximately twice the number of masers - not enough 
to account for the apparent deficit. 

Sites of maser emission observed with a single dish may be 
resolved into multiple sites when observed at higher angular res- 
olution with an interferometer. This will not affect the methanol 
masers counts as they all have interferometric follow-up, but may 
affect the HOPS water maser distribution. To investigate this, we 
compared the HOPS detections to ~4 x deeper, interferometric wa- 
ter maser observations t oward a 0.5 deg 2 region surrounding the 
GC dCaswelletai]|201ll) . In the same region (|/|, \b\ < 0.4°), 27 
masers were found compared to the 8 found by HOPS, and 3 of the 
8 HOPS detections were resolved into 2 maser sites. If this is repre- 
sentative of the whole GC-only region HOPS may have underesti- 
mated the number of water masers by a factor of 3. Preliminary re- 
sults from a deeper, interferometric EVLA water maser survey cov- 
ering \l\ < 1.5°, \b\ < 1° also recover a similar factor of 3 increase 
in number of detections compared to HOPS over the same region 
(J. Ott, F. Yusef-Zadeh private communication). This is not suffi- 
cient to explain the large disparity between the maser and dense 
gas distributions in the GC-only and non-GC regions. Just as im- 
portantly, there is also no reason to expect a systematic difference 



between the number of maser sites each HOPS detection will be re- 
solved into between the GC and non-GC regions. This is confirmed 
by preliminary results from the interferometric follow-up observa- 
tions of the HOPS masers (Walsh et al, in prep). As a further check, 
repeating the analysis using the integrated intensity of the masers 
(which will not be affected by how many maser sites each HOPS 
detection is resolved into at higher resolution), rather than number 
of masers, as a function of longitude produces the same results in 
Figures |2]& HJ. We adopt a conservative factor of 3 uncertainty in 
the maser counts, which we illustrate by the error bar in the middle 
panel of Figure|3] The similarity of both maser distributions in the 
bottom panel of Figure [3] gives confidence in the robustness of the 
maser counts and the adopted uncertainty. 

Finally, the trends in Figure [3] consist of ratios between ob- 
served parameters rather than the absolute values of the parameters 
themselves. The distance-dependence on expected source flux is 
the same for all observed parameters, so the trend itself will not be 
affected by systematic differences in distance between the GC-only 
and non-GC regions. We conclude that there is no evidence that our 
results may be strongly affected by observational biases. 

5.2 Potential systematic biases in the dense gas and SF 
tracers 

Many mechanisms can affect assumption i) above, the most im- 
portant of these being: the emission optical depth (assumed to be 
optically- thin); ammonia abundance and gas excitation conditions; 
dust properties (composition, temperature) and gas-to-dust ratio. 
These properties will certainly vary from one molecular cloud to 
another and undoubtedly cause much of the scatter seen in Fig- 
ure [3] as well as the non-linear slope in the NH3(1,1) integrated 
intensity vs 500/im flux. The variation in these properties must 
be understood before interpreting differences between individual 
clouds or deriving physical properties of the gas directly from one 
of these tracers. However, we are averaging over many molecular 
cloud complexes and are interested in order-of-magnitude varia- 
tions across the Galaxy so focus on systematic differences between 
GC-only and non-GC regions. 



5.2.1 Potential systematic biases in the dense gas mass estimates 

Regarding the mass estimates, the emission optical depth will on 
average likely be much higher towards the GC-only regions, lead- 
ing to a substantial underestimate of the dense gas mass towards 
the GC. This only accentuates the difference between the Galac- 
tic centre and the rest of the Galaxy. However, many observations 
show the GC excitation conditions and dust temperature are on 
average higher. The most relevant temperature measurements for 
this work are those derived directly from NH3. Targeted obser- 
vations of multiple NH3 inversion transitions towards GC GMCs 
show these have at least two kinetic temperature regimes - a cool 
(25 K) component which dominates the column density and a warm 
(200 K) component which acco unts for ~25% of the column den- 
sity dHuettemeister et al]| 19931) . Beam-averaged temperature mea- 
surements at the HOPS resolution results in an average Tk of 
~50K for GC clouds (J. Ott, priv comm.). Compared to ~20K 
for typical GMCs, this would lead to an overestimate of the GC 
dense gas mass by a factor 2 — 3. 

As the average gas kinetic temperature is higher in the GC, the 
NHs(l,l) effective critical density will be lower. Given the uniform 
sensitivity of the HOPS observations, this means the GC NHs(l,l) 
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integrated intensity will include emission from lower density gas. 
As the ultimate goal is to compare the amount of gas across the 
Galaxy above a single density threshold, including lower density 
material will lead to a systematic overestimate of this value towards 
the GC. The magnitude of the overestimate will depend on both the 
difference between the effective critical density probed towards the 
GC and the rest of the Galaxy, and how much additional gas in 
the GC lies between these two effective critical density limits. The 
NHa(l,l) effective critical density changes by less than a f actor 2 
for ga s with kinetic temperatures between 10 K and 100 K dEvansI 
1 19991) . The line brightness temperature of 1 K used in these cal- 
culations corresponds to 5 times the HOPS r.m.s. sensitivity so is a 
sensible limit. As the average gas kinetic temperature difference be- 
tween the GC and the Galaxy is much less than this (see above), the 
effective critical density probed will be very similar. In §[6] we esti- 
mate the average density of gas in the CMZ to be ~ 5 x 10 3 cm -3 , 
so most of the gas will be close to or above the effective critical 
density. We conclude systematic differences in the effective critical 
density probed between the GC-only and non-GC regions will not 
affect the result. 

Extrapolating the obser ved metallicity gradient in the disk of 
—0.03 to — 0.07dexkpc _1 (Balseret al. 201l]), one would expect 
the metallicity to increase by a factor 3—4 from the sun to the 
Galactic centre. While some studies do find an increased metal- 
licit y towards the Galactic centre, others find c lose to solar val- 
ues ^Shields & Ferlandl 19941 : iNaiarro et ail 2009). Using the com- 
monly adop ted assumption tha t the Galactic centre metallicity is 
twice solar (Ferr iere et alj2007h . the gas-to-dust ratio in the Galac- 
tic centre is likely lower by the same value. This means the 500/im 
flux towards the Galactic centre systematically over-estimates the 
dense gas mass by a factor ~2 comp ared to the disk. Simi larly, the 
NH3 abundance gradient in the disk (jDun ham e t al. 201 1) implies 
an average factor ~2 relative overestimate in the Galactic centre 
NH3-derived mass estimate. Both the 500/im and NH3 dense gas 
mass estimates therefore over-predict the Galactic centre mass by a 
factor ~2. No systematic effects are known regarding dust compo- 
sition with Galacto-centric radius. 

Another potential worry is related to beam dilution. The 
NH3(1,1) emission from the GC covers a large angular area on 
the sky and thus fills the 2' beam of the HOPS observations. If 
the molecular clouds outside of the GC systematically had angu- 
lar sizes much smaller than the 2' beam, the measured NHs(l,l) 
surface brightness would be affected by beam dilution and the inte- 
grated intensity would be systematically underestimated. To assess 
whether this affect is important, we estimated the distance at which 
molecular clouds of a given mass will be the same angular size as 
the 2' HOPS beam . The physical radius as a f unction of mass was 
estimated using the lKauffmann & Pillail J20101) empirical mass-size 
relationship, M = 87OM0 (r/pc) 1 ' 33 , for molecular clouds which 
will proceed to form high mass stars. Based on this relation, the line 
in the top panel of Figure[6] shows the distance at which a molecu- 
lar cloud of a given mass will have an angular size of 2'. Molecular 
clouds below the line will suffer from beam dilution so the observed 
surface brightness will be reduced, and the integrated intensity un- 
derestimated. From this we conclude that beam dilution is undoubt- 
edly an important factor for low-mass clouds at large distances. An 
implicit assumption in the above analysis is that it is possible to de- 
tect the emission from all the gas in the cloud. However, the density 
fluctuations in molecular clouds mean that some of the gas may lie 
at low density and not excite NH3 emission. We can estimate the 
effect of this by calculating the averag e density of the gas inferre d 
from the empirical mass-size relation teauffmann & Pillail koiO), 



n(H 2 ) = 3.13 x 10 3 (r/pc)' 1 67 cm" 3 . As shown in the bot- 
tom panel of Figure [6] clouds up to ~1OOOM0 have average vol- 
ume densities several 10 3 cm~ 3 and above, so we would expected 
most of the gas to emit in NH3. The average density of more mas- 
sive clouds drops, so the fraction of NH3 emitting gas will drop 
too. Although the top panel of Figure [6] shows HOPS will easily 
resolve very massive clouds to large distances, the NH3 emission 
must come from denser clumps within these clouds. The extent to 
which beam dilution affects the detected NH3 emission, if at all, de- 
pends on the density fluctuations within the clouds, which will vary 
from cloud to cloud. However, the fact that in general the larger the 
cloud, the better it will be resolved, mitigates the fact that less of the 
gas will be at high density. We conclude that while beam dilution 
undoubtedly affects low-mass clouds at large distances, it is not an 
important issue for clouds of a thousand solar masses and greater. It 
is the clouds in this mass regime that will be forming the high-mass 
stars traced by methanol maser emission and HII regions. 

On balance, the systematic effects would argue that making 
assumption i) would overestimate the dense gas mass towards the 
GC relative to the rest of the Galaxy. However, it seems highly un- 
likely these would lead to an overestimate by a factor of ~100. 
Additionally, as most of these effects are (at least partially) inde- 
pendent between the 500 /im and ammonia emission, the fact that 
no jump is seen in the correlation between the two tracers in Fig- 
ure[3]when comparing the GC-only and non-GC regions argues the 
systematics are not important for the analysis presented here. Based 
on the above arguments and the observed scatter in the 500 /im vs 
ammonia emission correlation (top panel of Figure [3}, we adopt a 
factor 5 uncertainty, which is illustrated by the error bars in Fig- 
ure [3] 

5.2.2 Potential systematic biases in the maser counts 

As water masers are also observed towards evolved stars, contam- 
ination may affect the maser number counts. As previously stated, 
the water masers from evolved stars are found outside the disk- 
plane and the emissio n is typically two orders of magnitude weaker 
jCaswelletal]|201(]|) . Therefore, sensitivity-limited, water maser 
Galactic plane surveys should predominantly trace SF regions. 
Nevertheless, we are conducting high angular resolution follow-up 
observations of the HOPS masers to uncover the true nature of each 
maser detection. Preliminary results show the majority (~70%) are 
associated with SF. Also, no trends in the relative number of masers 
associated with evolved stars or SF regions are seen with Galactic 
longitude (Walsh et al, in prep). In any case, if anything, one might 
expect more evolved stars towards the bulge and Galactic centre. 
Contamination of water masers from evolved stars would lead to 
an increase in number of water masers observed towards the GC. 
This only strengthens the result. Finally, the similarity of the water 
maser distribution in the bottom panel of Figure[3]to the methanol 
maser distribution (which are not observed towards evolved stars), 
adds further weight that contamination does not affect the results. 

Given the tight constraints on the physical conditions re- 
quired for masers to exist, another possibility is that maser emis- 
sion is somehow quenched towards the Galactic centre region. This 
seems plausible given the extreme environmental conditions in the 
Galactic centre compared to the rest of the disk (e.g. interstel- 
lar radiation field enhanced by 10 3 , external pressure P/k^lO 8 , 
strong magnetic fields, turbulence, gamma rays : lMorris & Serabvnl 
1 19961 : iFerriere et alj|2007l ; ICrocker et alj|2010h . However, masers 
are formed deeply embedded within GMCs so they are shielded 
from the interstellar radiation field. Regarding the effect of turbu- 
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lence, the coherence path lengths for masers are of order AU size 
scales, so individual maser emission regions will not be affected by 
velocity flows or gas kinematics on pc-scales (Mark Reid, private 
communication). Additionally, extremely bright "mega masers" are 
known to exist in the circumnuclear disks around supermassive 
black holes - an even more extreme environment than the CMZ. It 
is worth pointing out that the pumping mechanisms for water and 
methanol masers are different. As far as we are aware, no known 
effect in maser pumping theory can account for the observed fac- 
tor 100 deficit in the number of both water and methanol masers 
towards the GC-only region. We conclude there is no evidence that 
our results are strongly affected by systematic biases in the dense 
gas and SF tracers. 



6 QUANTITATIVE ANALYSIS OF THE SFR AND GAS 
MASS IN THE CMZ 

The previous qualitative analysis strongly points to large differ- 
ences between the SFR per unit mass of dense gas in the CMZ 
and the rest of the Galaxy. However, the inability to derive absolute 
SFRs from first principle for masers and the large uncertainties in 
deriving masses from NHa(l,l) integrated intensities preclude us 
from a more quantitative analysis using these tracers. Therefore, 
we now focus on the CMZ, the region which stands out as poten- 
tially discrepant, and aim to directly test the predictions of different 
SF relations. 



6.1 Structure and physical properties of the molecular gas in 
the Galactic centre 

The Galactic centre environment has be en well studied 
(Morris & Seraby^ 1 19961 ; iFerriere et al] [2007). Here we sum- 
marise the structure and physical properties of the molecular gas 
within a Galacto-centric radius (roc) of 500 pc (|Z| < 3.5°) 
- ie the central molecular zone (CMZ). The CMZ can be 
decomposed into an outer and inner component, separated at 
roc ~ 230 pc, and further separated the inner CMZ into a 
"disk" and torus, with approxima te radial ranges of t o 120pc 
and 130 to 230 pc, respectively dLaunhardt et alj |2002|). Base d 
on the far-IR dust continuum emission dLaunhardt et al] |2002l) , 
the mass of the disk and torus components are derived to be 
~4x 10 6 M Q and ~1.6x 10 7 M Q , respectively, with a total CMZ 
mass (tgc ^ 500pc) of ~6x10 7 Mq. Despite the many sys- 
tematic uncertainties in dust-derived masses, and the notoriously 
uncertain X-factor affecting CO-derived mass estimates, there is 
general agreement in the literature that the total molecular ga s 
mass within r G c < 500 pc is 2- 6x 10 7 M Q dFerriere et"alf 2007). 

Numerous studies have argued for two distinct molecular gas 
components: a warmer, low-density component (^10 3 cm -3 ) of 
"diffuse" clouds with a large filling factor and comprising ~30% of 
the mass, and a high-density (^10 4 cm~ 3 ) com ponent with a vol- 



ume filling factor of roughly a few percent (see lFerriere et alj2007l . 
for details). Given the effective critical density of the NH3(1,1), it 
is likely the HOPS observations are tracing the dense component. 

The recent far-IR HiGAL survey provides a new, high an- 
gular resolution vie w of the CMZ molecular gas properties 
dMolinari et alj|201lh . We derive the column density maps for the 
region —2.5° < I < 3.5°, \b\ ^ 0.5° using well-tested methods 
(Batte rsbv et alj|201ll) . A major systematic uncertainty is the as- 
sumed dust opacity and we note that we systematically underesti- 
mates the mass by a factor 2—3 compared to the DUSTEM method 



(e.g. lMolinarietalJl201lh for this reason. Figure [4] shows the to- 
tal mass of gas which lies above a range of different column den- 
sity thresholds. The different curves show the effect on the derived 
column density of using different wavelength bands and different 
source-extraction algorithms. From this we conclude the mass de- 
rived from the column densities to be robust to the background and 
source extraction to the 10-20% level. 

The vertical lines in Figure |4] show the lLada et al.1 d2012h 
proposed column density th reshold of 7.5 x 10 21 cm -2 , calculated 
from the lLada et al] d2012l ) extinction threshold of Ay = 8 mag 
and assumi ng an Ay — » Nh-, co nversion of Nh 2 = Av x 0.95 x 
10 21 cm -2 dFrerking et al.1 1982b . This shows that most of the mass 
in this region lies above an extinction of Av = 8 mag. The masses 
derived from this analysis are reported in Table [2] The total de- 
rived mass of the CMZ gas from the HiGAL data is 4.1 x 10 7 M , 
in good agreement with previous measurements in the literature 
dFerriere et al.ll2007h . 

Determining the volume density of the gas from the HiGAL 
column density maps is not straightforward as it depends on knowl- 
edge of the 3D gas structure. For the most part this is not well con- 
strained towards the Galactic centre. The region of the CMZ with 
the be st quantified 3D stru cture is the molecular torus, or "100 pc 
ring" dMolinari et al]|20Tl|) . The volume of gas in a ring of radius 
R and width AR, with a thickness of H, is 



V = 2nRARH = 4.4x10 



60 / R 



VlOOpc 



AR 



12 pc J V20pc 



H 



(1) 

The fiducial radius of lOOpc is that determined for the ring and 
the fiducia l width of 12pc is b ased on the observed column den- 
sity maps (Molinari et al. 2011). We have used H — 20 pc assum- 
ing hydrostatic equilibrium (H/R — Vt/v c )', the circular velocity 
v c ~ lOO kms -1 while the turbulent velocity vt = \/8 In 2a v — 
20kms _1 dOka et alj200lh . 

With a mass of 1.8xlO 7 M (see column 5 of Table |2j, the 
volume density for the ring is then 



rripV 



5000 



( R 



VlOOpc 



AR 
12 pc 



H 
20 pc 



(2) 

Not e that this mass is a factor ~ 1 .7 smaller than reported pre- 
viously (Moli nari et al.l 120 111) due to the different dust opacities 
used in that work. The inferred densities in this work are there- 
fore also lower by the same factor. The surface density of the ring 
is 



M a 



2-kRAR 



2400 



R 



AR 



100 pc ) \12pc 



M 



©PC 



0.5 gem 



(3) 



(4) 



As the gas in the ring is observed to be clumpy, the density at in- 
dividual positions along the ring may clearly deviate from these 
representative values. 

The gas in the ring has been modelled as a twisted ellipse 
with s emi-major and semi-m inor axes of 100 pc and 60 pc, respec- 
tively (Mo linari et alj|20lll) . The dashed vertical lines in Figure [5] 
show the average volume density assuming the ellipse has the av- 
erage cylindrical diameter (ie H = AR) specified by the anno- 
tated values. From the measured height of the ring projected on 
the plane of the sky we estimate the average volume density to be 
~ 5 x 10 3 cm -3 . This agrees well with the analytical estimate 
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above and also volume densities estimated by previous authors for 
the dense component of the CMZ molecular gas (see above). 

We note that the gas volume density and interpretation of the 
measured column density changes considerably if the gas structure 
is different than assumed above. For example, in the worst-case 
scenario, if the gas assumed to make up the ring actually lay in a 
uniform disk with the same radius, the volume density would drop 
to around a few 10 2 cm" 3 . Additionally, the measured column den- 
sity would be highly dependent on the viewing angle, and appear 
much larger if viewed through the disk-plane rather than in a face- 
on view of the Galaxy. In this case it would not be fair to test the 
predictions of the column density threshold SF relations using the 
directly measured column density without taking this into account. 

However, several lines of evidence suggest this is not the 
case. On top of the mo delling which succe ssfully reproduces the 
gas kinematic structure (Molina ri et alj2oTlh . analysis of emission 
from many high critical density molecular line transitions show the 
majority of the gas in the CMZ is cold, has a density >10 4 cm" 3 
and a volume filling factor of only a few percent. In addition, if 
the gas was uniformly distributed, one would expect optically thin 
emission (e.g. from sub-mm dust emission) to peak along the line 
of sight with the largest path length, i.e. directly through the centre. 
In fact the emission intensity increases towards the edges, directly 
contradicting what would be expected from a uniform disk, and as 
expected from the emission arising in a ring. 



6.2 Measured CMZ SFR values 

Values of the SFR for the CMZ in the literature vary by fac- 
tors of several, but recent measurements based on number counts 
of massive young stellar objects, appear to converge to a value 



close to 0.1 M, 



lYusef-Zadeh et all2009l ; llmmer et al.ll2012l) . 

As an alternative approach, we make use of the recent anal- 
ysis of WMAP data dMurrav & Rahmarj|20ld) to determine the 
fraction of ionising photons in the Galaxy coming from the GC. 
This analysis offers many advantages over older estimates of the 
ionising photon flux. WMAP has much better absolute flux pre- 
cision, suffers less from spatial filtering and observed at multi- 
ple frequencies so the spectral index of the emission can be cal- 
culated allowing the relative contributions from free-free, non- 
thermal and spinning dust emission to be determined. Although the 
WMAP beam is large, ~ 1° at the lowest frequencies, the free-free 
emission in the Galaxy is dominated by the so-called "extended 
low density" (ELD) Hn e mission, rather than classical Hn regions 
dMurrav & Rahmarj|20 1 fj) . 

Based on the free-free-only contribution to the 33 GHz flux, 
corrected for dust attenuation, the rate of ionising ph otons, Q, 
acros s the Galaxy is measured to b e 2.9xl0 53 s~ 1 iLeeetall 
l2012h . The SFR, M», is then given by dMurrav & Rahmarj201Ch 



(e.g. O.MMflyr" 1 and 0.08 M yr- 



the IMrfl Murray & RahmarJ d2010h show that 

(m„) 



(?) 
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(t > =3.9 x 10 6 yr- 



(6) 



(7) 



M. = Q 



ii) (t Q y 



(5) 



The reason this timescale is small compared to the average lifetime 
of stars in a cl uster is because {to) is dominated by stars >40M© 
dMurray||201lh . While stars of ^8 Mq (with lifetimes significantly 
longer than (tq)) produce ionising photons, the main-sequence 
life time of stars >40Mq is a slowly varying function of mass, 
while the ionising photon output continues to increase rapidly. The 
WMAP observations are therefore sensitive to star formation over 
the last ~4 Myr. 

From the total rate of Lyman continuum photons in the 
Galaxy of Q =2.9x 10 53 s" 1 , the tot al SFR is calculated to be 
1.3 Moj/r" 1 dMurrav & RahmarfeOlOl) . In Table|2]we show the Q 
values for the same latitude and longitude regions as the previous 
gas mass calculations. The range of values in column 7 show the 
total Q assuming the near and far kinematic distances. Column 9 
gives the fraction of the total Galactic Q that the average of these 
values corresponds to. Based on this analysis, for the latitude range 
of |6| < 0.5°, the longitude ranges -2.5° < I < -1.0°, \l\ < 1° 
and 1° < I < 3.5° therefore have SFRs of ~0.016, ~0.015 and 
-~0.0036Mq yr" 1 , respectively. We report these numbers in Ta- 
ble |2] 

This means that in total ~2.8% of the star formation in the 
Galaxy lie in this / and b range. Although this encompasses the 
region most people would use to define the CMZ, given the differ- 
ent angular resolutions of the studies used to derive the mass and 
SFR, we check to make sure the hard I and b boundary imposed 
is not affecting this result. The WMAP data extend to higher lat- 
itudes, and we note that including the latitude range \b\ < 1° for 
the \l\ < 1° region increases the maximum Q by a factor ~3. The 
total SFR over this larger region is 0.06 Mg yr~ . The Q values 
over the larger latitude range are shown in parentheses in columns 
7, 9 and 10 of Table[2] The other longitude regions are not affected 
by increasing the latitude range. As a final check, we note that the 
total Q in the much larger region |/| < 4°, |6| < 1° is 1.64x 10 52 , 
corresponding to 5.6% of the total Galactic SFR. 



6.3 Comparing measured CMZ gas mass and SFRs with 
those predicted from SF relations 

We now compare the measured gas mass and SFRs determined for 
the CMZ above, with the values predicted from the column-density 
threshold and volumetric SF relations. 



6.3. 1 Column-density threshold SF relation: predicted CMZ SFR 

The lLada etafl J2012h column-density threshold relation proposes 
the SFR in a molecular cloud is set by the amount of gas above 
a column density threshold corresponding to an extinction of 
Av = 8 mag. This is parametrised through, SFR = 4.6 x 



where (q) is the ionising flux per star averaged over the stellar ini- 
tial mass function and (m*) is the mean mass per star. The quantity 
(to) is the ionisation- weighted stellar lifetime, or time at which the 
ionising flux of a star falls to half its maximum value, averaged over 



3 Assuming high mass stars are well selected from the IMF as expected if 
their formati on is intrinsically linked with the formation of massive stellar 
clusters (e.g. Smith et al . 2009). However, this assump tion may be affecte d 
if stochastic sampling of the IMF is important (see e.g. Br essert et alj2012h . 
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10~ 8 (/d) Mgas(M Q )M Q yr~\ where the factor f D is the frac- 
tion of the total gas mass (M gas ) above the column density thresh- 
old. Table|2]shows the mass of gas above the column density thresh- 
old for the three different CMZ regions. Based on these mass mea- 
surements, the predicted SFRs for these regions are 0.18, 0.78 and 
0.74 M©yr -1 . These values are listed in column 1 1 of Table[2] 

6.3.2 Volumetric SF relation: predicted CMZ SFR 

Next we consider the Rrumholz et alj d2012l) volumetric star for- 
mation relation which is parametrised through: 

SFR M mol 1 M mol 1 

[volume] [volume] rg [volume] Ts 

where the SFR in a given volume of gas per free-fall time (rs) is 
determined by the mass of the molecular gas (M mo i) in that volume 
and a global efficiency of gas converted to stars per free-fall time, 
rj, of 1%. The free-fall time is 

The free-fall time, and hence gas density, for the gas in the 
GC clearly plays an important role in setting the predicted SFR. As 
discussed in § 16.11 the region of the CMZ with the best constrained 
gas structure is the 100 pc ring (see §[67TJ. From Eq[8]&|9] the SFR 
predicted by the volumetric relation is given by, 

M*=vM a /r ff = .53( m ^) M yr"\ (10) 

for r\ = 0.01. This is twenty times the observed star formation 
rate. The fi ducial predicted SFR increases to 0.9 Mq yr _1 for rj = 
0.017 (e.g. lKennicutll998h . 

Given the clumpy nature of the gas in the ring, for illustra- 
tive purposes, Figure[5]plots the predicted SFR as a function of the 
gas density. The dash ed lines show the range in observed thickness 
(Moli nari et al.l20l"u) . or scale height H (and hence the inferred ra- 
dial extent AR if the ring is cylindrical) of the ring from ~6pc to 
factors of several larger. 

We note that, if instead of using the observed geometry of gas, 
one assumes the GC gas mass is spread over a uniform disk of di- 
ameter ~100pc, the expected free-fall time significantly increases, 
and the predicted star formation rate drops considerably. Using this 
much larger free-fall time might be appropriate if the timescale for 
star formation in the GC was set by the time it would take the 
gas to collapse to the centre of the ring. However, the observed 
gas dynamics show that the ring is not collapsing to wards the cen- 
tre, b ut rather rotating around the centre of the ring (Molin ari et al.1 
Urn]). In addition the gas in the ring clearly fragments, and the 
free-fall time at the observed local density is much shorter that the 
time it would take to collapse to the centre (assuming no rotational 
support). We conclude that the observed local volume density of 
~ 5 - 8 x 10 3 cm -3 is the correct density to use for the purpose 
of calculating the free-fall time. 

We believe the reason our results contradict previous au- 
thors who have concluded the CMZ lies on the S -K relation (e.g. 
lYusef-Zadeh et alJl2009l : lKenmcutt & Evandl2012l) is because they 
assumed the observed CMZ gas lay in a uniform disk and thus had 
a much lower average volume density. 

Assuming the rest of the gas in the outer CMZ is at a similar 
density to that in the torus l lFerriere et al.ll2007h . the predicted SFR 
for the -3.5° < I < -1° and 1° < I < 2.5° regions are 0.14 and 



0.39 M©yr , respectively. The predicted SFR values are listed in 
column 12 of Table[2] 



6.3.3 Predicted vs measured SFR values 

The predicted SFRs in Table|2]for both the column density thresh- 
old and volumetric SF relations agree to within a factor of two. The 
measured SFRs are at least a factor 10 smaller than the predicted 
values for all three regions. The region of the outer CMZ at positive 
longitudes (ie 1° < I < 3.5°) stands out as extreme in this regard. 
Despite containing almost the same reservoir of gas as the inner 
CMZ (ie \l\ < 1°), it has a much lower measured SFR. Both the 
column density threshold and volumetric SF relations over-predict 
the SFR in this region by two orders of magnitude. 

However, the SFR predictions for the volumetric SF relation 
in Table [2] rely on the assumption that the volume density of the 
outer CMZ is similar to that determined for the inner CMZ "100 pc 
ring". Given the uncertainties in determining the gas density from 
the measured surface density, we can approach the problem from 
a different direction and ask, what gas densities would be required 
to produce the measured SFR from the measured gas masses? As 
shown in Table[2] the 1° < I < 3.5° region contains 0.85% of the 
total molecular gas in the Galaxy and 0.3% of the total SF. From 
Eq[8] the density of the clouds in this region should be 0. 1 times the 
density of those in the disk, or 0.1 x 250 cm~ 3 ~ 25 cm -3 . This is 
unfeasibly low. 

We conclude the column density and volumetric SF relations 
over-predict the SFR in the CMZ by an order of magnitude. 



7 SUMMARY AND DISCUSSION 

In summary, we find that the dense gas (NH3 & 500 /im) and SF 
activity tracers (masers & HII regions) used in this study are reli- 
ably tracing the present-day relative dense gas mass and SF activity 
distributions, respectively, across the Galaxy. We conclude that the 
striking difference between the dense gas and SF activity tracers 
between the GC-only and non-GC regions shows the current star 
formation rate per unit mass of dense gas is an order of magni- 
tude smaller in the GC than in the rest of the Galaxy. We directly 
test the predictions of proposed column-density threshold and vol- 
umetric star formation relations and find, given the mass of dense 
gas in the GC, they over-predict the observed SFR in the GC by 
an order of magnitude. We conclude the current star formation re- 
lations are incomplete in some way. Any universal column/volume 
density relations must be a necessary but not sufficient condition 
for SF to occur. 

Putting the CMZ in the context of the Galaxy as a whole, the 
Milky Way contains ~ 2 x 10 9 M of gas, so the CMZ holds 
roughly a few percent of this. The WMAP analysis shows the CMZ 
also contains a few percent of the ionising photons, and hence star 
formation, in the Galaxy. Both the volumetric and surface density 
SF relations predict that a given mass of gas will form stars more 
rapidly if the respective densities are larger. Yet the gas in the CMZ, 
which has a much higher surface and volume density than any com- 
parable mass of gas in the disk of the Milky Way, forms stars at a 
rate proportional to the ratio of gas in the CMZ to that in the Milky 
Way. 

Something is required to slow the rate of SF in the CMZ com- 
pared to that in the rest of the Milky Way. An additional support 
mechanism, not taken into account in either the column density 
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threshold or volumetric SF relations, may be responsible for in- 
hibiting the SF. One potential solution is therefore an additional 
term or threshold in the proposed SF relations. The most notice- 
able difference between clouds in the CMZ and the rest of the 
Milky Way, apart from the 1 to 2 order of magnitude larger vol- 
ume density, is the order of magnitude larger internal cloud veloc- 
ity dispersion ^Morris & Serabvn|[l996l ; iFerriere et al.ll2007h . We 
thus conclude it is likely that the relevant term is related to the 
additional turbulent energy in the gas providing support against 
gravitational collapse. A simple way to parametrise this is through 
the linewidth ratio, AK- a tio = AV/AVo, where AVo is the typ- 
ical internal cloud velocity dispersion in disk molecular clouds. 
The Schmidt law could then be re-expressed either in terms of 
surface density, T, S fr oc (E gas ) Q /(AKatio) i ' [a ~ 1 to 1.4; 
b ~ 1], or mass of gas above the density threshold, Md enS e, 
through, J^sfr oc Afd C nsc/(AV r atio) t '. This would reconcile the 
star formation in the CMZ with that in t he rest of the Galaxy. Theo- 
rists have wrestled with this before (e.g.lKru mholz & McKe el2005l : 
IPadoan & Nordlundll201lUKrumholz et alj|2012h and we are cur- 
rently seeking to test these scenarios. 

However, even if the extreme environmental conditions in the 
CMZ do inhibit SF, they do not stop it entirely. The CMZ con- 
tains Sgr B2, one of the most extreme cluster forming regions in 
the Galaxy. It also contains high-mass star clusters like the Arches 
and Quintuplet, and at least one molecular cloud which appears to 
be the progenitor of such massive stellar clusters dLongmore et al.l 
l2012l ;lBressert et al. 20 I2T). Evidence exists of episodic SF events in 



the CMZ dSofue & Hand J 19841; 



l20ld : iBlan d-Hawthorn & Cohen 



Yusef-Zadeh et alJl2009l ; ISu et al.l 
20031) and mechanisms exist to 



explain how such episodic SF can occur. Gas in barred spiral galax- 
ies like the Milk y Way is funnelled from the disk through the bar 
to the GC dKormendv & Kennicuttll2004l ; ISheth etafll2005t) . If gas 
is continually fed from the disk to the GC, and the environmental 
conditions impose a higher threshold for SF to occur, the gas might 
build up until reaching a critical point before undergoing a burst of 
SF. 

Although clouds near the Galactic centre in the Milky Way 
and other galaxies may only represent a small fractional volume 
of a galaxy, they can contribute a significant fraction of the total 
dense molecular gas. In terms of dense gas mass and environmental 
conditions, the Galactic centre also acts as a bridge between local 
SF regions in our Galaxy and SF environments in external galaxies. 
Understanding why such large reservoirs of dense gas deviate from 
commonly assumed SF relations is of fundamental importance and 
may help in the quest to understand SF in more extreme (dense) 
environments, like those found in interacting galaxies and at earlier 
epochs of the Universe. 
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Figure 1. Distribution of dense gas and star formation activity tracers as a function of Galactic latitude and longitude. Black circles mark the positions of 
Hn regions, blue crosses show methanol masers and red plus symbols mark water masers. Regions of sky not covered in these surveys are shaded in grey. 
NH3(1,1) [bottom] and H69a [second-bottom] integrated intensity emission is displayed using a square-root image stretch. The tracers and their function as 
either a dense gas or star formation activity tracer are labelled at the right hand edge of the bottom row. The central molecular zone (CMZ) can be seen as 
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Figure 2. [Top] Dense gas distribution, traced by the NH3(1,1) integrated intensity, as a function of Galactic longitude (2' bins). The emission was split into 
two regions: a Galactic centre region encompassing the central molecular zone (blue) and the rest of the Galaxy (green). The Galactic centre region clearly 
dominates, accounting for ~80% of the NH3(1,1) integrated intensity emission, despite only accounting for ~4% of the total survey area. The inset shows 
the same plot but with a logarithmic scale on the y-axis. [Bottom] Star formation activity, traced by the number of water masers (grey) and methanol masers 
(black: data for I >20° unavailable), as a function of Galactic longitude (0.25° bins). In comparison to the NH3(1,1) emission, the distribution of both maser 
species is relatively flat with Galactic longitude. 
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Figure 3. [Top] Comparison of two independent dense gas tracers, NH3(1,1) integrated intensity vs. the integrated 500 /im flux, for each 2° longitude bin 
(normalised to the total flux in that tracer). [Middle] The fraction of dense gas tracers (NH3(1,1) integrated intensity and 500/^m flux) vs fraction of star 
formation tracers (water and methanol masers) for each 2° longitude bin between —70° < I < 20° (the mutually-overlapping coverage of MMB and HOPS). 
Different colours show different combinations of dense gas and SF tracers. In both top and middle plots, the Galactic centre (GC-only) and disk (non-GC) 
regions are shown as triangles and crosses, respectively. Representative uncertainties of a factor 5 in dense gas mass and 3 in maser counts are shown as error 
bars. [Bottom] Ratio of dense gas surface density (NH3(1,1) integrated intensity) to star formation surface density (methanol and water masers) as a function 
of Galacto-centric radius, Rgc- The hatched region shows the radial range over which the rotation curve used to derive Rgc i s unreliable. The surface density 
ratios over this region should be ignored. 
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Figure 4. Total mass in various longitude ranges of the CMZ which lies above a range of column density thresholds. The vertical dashed line shows the 
proposed column density threshold of 7.5 X 10 cm , calculated from the extinction threshold of Ay = 8 mag and assuming an Ay — > Nh 2 conversion 
of N[i 2 = Ay X 0.95 X 10 21 cm -2 . The different curves show the effect on the derived column density of using different wavelength bands and different 
source-extraction algorithms. From this we conclude the column densities are robust to the 10-20% level. 
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Table 1. Properties of the Galactic plane surveys used in this work. 



Survey 


Angular 
Resolution 


Transition or 
Continuum 


Frequency or 
wavelength 
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Notes 
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(s„ km^;") 



Table 2. Gas mass and ionising photon rate (Q) for various longitude and latitude ranges of the CMZ. The second row corresponds roughly to the longitude 
range for the inner ~150pc of the CMZ. The first and third rows correspond roughly to the longitude ranges of the outer ~ 150 to ~500pc of the CMZ, in the 
negative and positive longitude ranges, respectively. The gas masses are derived from the HiGAL column density maps. The value in column 5 is the total mass 
in that longitude and latitude range, whereas column 6 reports the mass above the column density threshold corresponding to an extinction of Ay = 8 ma g 
(see text for details). The range in values in column 7 corresponds to the total Q assuming near and far kinematic distances, respectively I Lee et al. 2012). 
Column 8 gives the fraction (in percent) of the total Galactic molecular gas contain ed in column 5, a ssuming a total Galactic molecular gas mass of2xlO 9 M . 
Column 9 gives the fraction (in percent) of the total Galactic Q of 2.9 X 10 53 s" 1 iLee et al.ll2012l) . Column 10 gives the corresponding SFRs for each of the 
regions. The values in parentheses in columns 7, 9 and 10 are the Q values over the same longitude regions if the latitude range is doubled to |b| < 1°. Only 
the |Z| < 1° longitude range is affected. The SFRs in columns 11 and 12 are those predicted from the column density threshold and volumetric SF relations, 
respectively (see text for details). 
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Figure 5. The predicted star formation rate from the Krumholz et al. 12012) volumetric relation as a function of the average gas volume density. The dashed 
vertical lines show the average volume density assuming the 1.8 X 10 7 Mq of gas in the region \l\ < 1°, \b\ < 0.5° (see Table|2j lies in a cylindrical ring 
with semi-major and semi-minor axes of 100 pc and 60 pc, respectively, with the diameter specified by the annotation for each line. The observed thickness 
(and hence the inferred diameter) of the ring varies from ~6pc to factors of several larger. We estimate the average volume density to be ~ 5 X 10 3 cm -3 , 
with a corresponding predicted SFR of ~0.4 M0/yr. 
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Figure 6. [Top] Estimating the distance at which molecular clouds of a given mass will be the same angular size as the 2' HOPS beam. The physical radius as 
a function of mass was estimated using the lKauffmann & Pillail fcoid) empirical mass-size relationship for molecular clouds which will proceed to form high 
mass stars. The line in the plot shows the distance at which a molecular cloud of that mass will have an angular size of 2'. Molecular clouds below the line 
will suffer from beam dilution so the observed surface brightness, and therefore integrated intens ity, will be reduced. Molecu lar clouds above this line will not 
suffer from this affect. [Bottom] Density of clouds as a function of their mass, calculated from the Kauffmann & Pillai 12010) empirical mass-size relationship. 



